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Abstract. We describe how the study of resonances and fluctuations can help constrain the thermal and chem- 
ical freezeout properties of the fireball created in heavy ion collisions. This review is based on I1I2I3I4I5I 



1 Introduction 

The idea of modeling the abundance of hadrons using sta- 
tistical mechanics techniques has a long and distinguished 
history [6 7 8 9I10I . In a sense, any discussion of the ther- 
modynamic properties of hadronic matter (e.g. the exis- 
tence of a phase transition) requires that statistical mechan- 
ics be applicable to this system ( through not necessarily at 
the freeze-out stage). 

That such a model can describe quantitatively the yield 
of most particles, including multi-strange ones, has in fact 
been indicated by fits to average particle abundances at 
AGS,SPS and RHIC energies I10I1 1I12I13I14I15I16I . The 
jump from noting the qualitative goodness of this descrip- 
tion to ascertaining its limits, and using it as a tool to study 
dense hadronic matter, is however a still ongoing process. 

Some practitioners [6| interpret the statistical model 
results in terms of nothing more than phase space domi- 
nance: For a process strongly enough interacting with enough 
particles in the final state, dynamics "factors out" into a 
normalization constant, and the final state probabilities are 
dominated by phase space. If this is the case, the applica- 
bility of the statistical model has nothing to do with a gen- 
uine equilibration of the system. Others think that in soft 
QCD processes particles are "born in equilibrium" [16], 
and the applicability of the statistical model to even smaller 
systems is a fundamental characteristic of QCD. Still oth- 
ers 01 11171 believe that the applicability of the statistical 
model is a sign of a phase transition, as the chemical equili- 
bration of hadrons signals a regime in which multi -particle 
processes and high-lying resonances dominate. 

Obviously, the link between the QCD phase diagram, 
defined in the Grand Canonical limit, and experimental 
data can only be made in the last interpretation. Yet it is not 
clear how these interpretations can be differentiated. Cur- 
rently, the debate centers around the scope of application 
of the statistical model [[1611 81 to smaller systems. Yet in 
every one of these cases an objective criterion for linking 
the goodness of fit to the effective relevance of the under- 
lying theory is still absent fl6l : Since it is clear that,in all 
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regimes, non-statistical processes such as jet fragmenta- 
tion and "corona physics" are also present, and since error 
bars vary dramatically across system sizes for essentially 
experimental reasons, excluding the statistical model from 
a raw x 1 analysis is highly nontrivial II 1 61 1 81 . 

Two, related, observables which might lead to progress 
in this context are short-lived QCD resonances [2.3 1 and 
event-by-event fluctuations 141191 . 

Short-lived resonances go to the heart of this question 
because,in principle, their abundance can be modified by 
elastic interactions of already formed hadrons. The abun- 
dance of hadrons stable against strong decays, such as the 
A, is expected to be unchanged even if particles created at 
hadronization continue to interact, because inelastic inter- 
actions (such as An — > pK) are suppressed w.r.t. purely 
elastic scattering. Even purely elastic scattering, such as 
An — > E* — > An could in principle create new E*s or 
make £*s created at hadronization undetectable (since such 
particles are detectable by invariant mass reconstruction 
only). Since the fitted hadronization temperature of T = 
170 MeV generally implies a non-negligible reinteraction 
phase for A-A collisions, but not for smaller systems, the 
applicability of the statistical model to describe unstable 
resonances as a function of system size is a good gauge to 
see to what extent the statistical model applies specifically 
at the hadronization stage. 

As for fluctuations, it is a fundamental principle of sta- 
tistical mechanics that variances around averages scale w.r.t. 
averages in a way defined by the maximization of entropy 
under the constraints specific to the ensemble. In our con- 
text "Averages" are particle multiplicities per event and 
fluctuations are event-by-event fluctuations. For macroscopic 
systems, this principle ensures that fluctuations become 
negligible and the expectation that the state of the system is 
the maximum entropy one is nearly certain to be realized. 
(9(100- 1000) particles is not enough for this to be the 
case, but, if statistical mechanics applies, one should still 
see that yields, fluctuations and higher cumulants scale in 
a way calculable from the partition function. 

These two classes of observables are actually related 
1 5 20.21 j: Even invisible resonances in the medium, whose 
decay products are rescattered and thus undetectable by 
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invariant mass reconstruction still continue to "exist" as 
a correlation between their decay products abundances. In 
contrast, "regenerated" resonances do not change such cor- 
relations , because they are formed by particles already 
generated at hadronization. Comparing resonance abundances 
with fluctuation observables can be used to gauge, from 
experimental observables, whether hadronization is where 
particle abundances are fixed once and for all, or whether 
further dynamical evolution also impacts such abundances. 

In the statistical model there are two types of chemi- 
cal equilibrium 1 1 1 . : all models assume relative chemical 
equilibrium, but some also assume absolute chemical equi- 
librium, implying the presence of the "right" abundances 
of valance up, down, and strange quark pairs to be present 
after hadronization. In absolute chemical equilibrium at 
highest heavy ion reaction energy one obtains chemical 
freeze-out temperature T ~ 160- 170 MeV, which goes 
down to ~ 50 MeV at the lowest reaction energies irTzl . 

The energy dependence of the freeze-out temperature 
1 1 1 than follows the trend indicated in panel (a) of figure 
[T] as the collision energy increases, the freeze-out tem- 
perature increases and the baryonic density (here baryonic 
chemical potential jUb) decreases IfTSl . An increase of freeze- 
out temperature with -\fs is expected on general grounds, 
since with increasing reaction energy a greater fraction of 
the energy is carried by mesons created in the collision, 
rather than pre-existing baryons |Q]. 

Further refinements in the approach described above 
are often implemented, the more notable ones are allowance 
for strangeness chemical nonequilibrium ITTJI at low -\[s 
and canonical effects for small reaction volumes El . These 
effects do not materially alter the behavior of temperature 
and chemical potential shown in the panel (a) of Fig. Q] 
Such smooth variations fail to fully reproduce the non- 
continuous features in the energy dependence of hadronic 
observables, such as the "kink" in the multiplicity per num- 
ber of participants and the "horn" [12.23.24 1 in certain par- 
ticle yield ratios. 

Non-monotonic behavior of particle yield ratios could 
indicate a novel reaction mechanism, e.g. onset of the de- 
confinement phase [ 23 1 . In a rapidly expanding fireball un- 
dergoing a change in microscopic degrees of freedom, it is 
a possibility that absolute chemical equilibrium does not 
hold. Either super-cooling [ 10] or bulk viscosity-triggered 
instabilities [25 26] could justify such lack of chemical 
equilibrium on microscopic grounds. 

One can model this phenomenologically by removing 
the hypothesis of absolute chemical equilibrium among hadrons 
produced. The systematic behavior of T with energy in this 
case is quite different |24l , as is shown in panel (b) of fig- 
ure Q] The two higher T values at right are for 20 (lowest 
SPS) and (most to right) 11.6 A GeV (highest AGS) re- 
actions. In these two cases the source of particles is a hot 
chemically under-saturated (T ~ 170 MeV ) fireball. Such 
a system could be a conventional hadron gas fireball that 
had not the time to chemically equilibrate. 

At higher heavy ion reaction energies it is possible [ 10] 
to match the entropy of the emerging hadrons with that 
of a system of nearly massless partons when one consid- 
ers supercooling to T ~ 140 MeV (essentially bringing 
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Fig. 1. (Color online) Der^ftiencO of freeze-out temperature T 
and baryo-chemical potential fi% on reaction energy in the Equi- 
librium (panel (a), 11121 ) and non-equilibrium (panel (b),[ 24 1) 
freeze-out. The arrow corresponds to increasing ^fs. 



chemical and thermal freeze-out close together, as seen in 
1 27 28 1), while both light and strange quark phase space 
in the hadron stage acquire significant over-saturation with 
the phase space occupancy y q =u,d > 1 an d at higher en- 
ergy also y s > 1. A drastic change in the non-equilibrium 
condition occurs near 30 A GeV, corresponding to the dip 
point on right in panel (b) of the figure Q] (marked by an 
asterisk). At heavy ion reaction energy below (i.e. to right 
in panel (b) of figure [TJ of this point, hadrons have not 
reached chemical equilibrium, while at this point, as well 
as, at heavy ion reaction energy above (i.e. at and to left 
in panel (b) of figure [T), hadrons emerge from a much 
denser and chemically more saturated system, as would 
be expected were QGP formed at and above 30 A GeV. 
This is also the heavy ion reaction energy corresponding to 
the "kink", which tracks the QGP's entropy density (higher 
w.r.t. a hadron gas), and the peak of the "horn" 11231 . which 
tracks the strangeness over entropy ratio (also higher w.r. t. 
a hadron gas). 

The main reason for the wider acceptance of the equi- 
librium approach y, = 1 is its greater simplicity, there 
are fewer parameters. Moreover, considering the quality of 
the data the non-equilibrium parameter y q is not necessary 
to pull the statistical significance above it's generally ac- 
cepted minimal value of 5 %. On the other hand, the pa- 
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rameters y q and y s were introduced on physical grounds 
1 10 1, thus these are not arbitrary fit parameters. Moreover, 
these parameters, when used in a statistical hadronization 
fit, converge to theoretically motivated values. They also 
help to explain the trends observed in the energy depen- 
dence of hadronic observables. 



2 Resonances 

Many strong interaction resonances, a set we denote by the 
collective symbol N* (such as K*°(892), ,4(1232), £*( 1 385), 
A*(1520), 3* (1530)) carry the same valance quark con- 
tent as their ground-state counter-parts N (corresponding: 
K, p, 27, A, 3). N* typically decay by emission of a pion, 
N* — > N +n. Considering the particle yield ratio (N*) / (N) 
in the Boltzmann approximation (appropriate for the par- 
ticles considered), we see that all chemical conditions and 
parameters (equilibrium and non-equilibrium) cancel out, 
and the ratio of yields between the resonance and it's ground 
state is a function of the masses, and the freeze-out temper- 
ature, with second order effects coming from the cascading 
decays of other, more massive resonances 0101151 : 



Tn) 



(1) 



where W(x) = x 2 K2(x) is the (relativistic) reduced one par- 
ticle phase space, K2(x) being a Bessel function, g is the 
quantum degeneracy, and b y« is the branching ratio of res- 
onance j decaying into R. 

Because of the radically different energy dependence 
of freeze-out temperature in the scenarios of [ 12 and [24], 
seen in figure [1] the prediction for the resonance ratios Eq. 
(UJ vary greatly between these two scenarios. In the equi- 
librium scenario the temperature goes up with heavy ion 
reaction energy, and thus the resonance abundance should 
go smoothly up for all resonances. On the other hand, the 
nonequilibrium scenario, with a low temperature arising 
only in some limited reaction energy domain, will lead to 
resonance abundance which should have a clear dip at that 
point, but otherwise remain relatively large. 

We have evaluated several resonance relative ratios shown 
in figure [2] within the two scenarios, using the statistical 
hadronization code SHARE II15I19H . For the non-equilibrium 
scenario, we have used the parameters given in [24], table 
I. For the equilibrium scenario, we used the parametriza- 
tion given in [12] figures 3 and 4. In the latter case, the 
strangeness and isospin chemical potentials were obtained 
by requiring net strangeness to be zero, and net charge per 
baryon to be the same as in the colliding system. 

As seen in figure [2] the expected trend with -\fs is ap- 
parent in all considered resonance ratios, though in cases 
where the mass difference is large, the effect is much more 
pronounced than in some others. Indeed, for many of the 
ratios we present the experimental error may limit the use- 
fulness of our results However, because of the cancellation 
(to a good approximation) of the baryo- and strangeness 
chemical potentials, the qualitative prediction for the en- 
ergy dependence of the resonance yields within the two 
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Fig. 2. (Color online) Ratio of resonance to the stable particle. 
Thick lines for particles with strange quark content, thin lines for 
particles with anti-strange quark content, as a function of energy. 
Solid black lines refer to the equilibrium fits (y q s = 1), with the 
parameters for AGS and SPS energies taken from 0121 . Dashed 
red lines refer to non-equilibrium fits (y q s fitted), with the best fit 
parameters for AGS and SPS energies taken from 0241 . 
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and the (N*) / (AO yield from that point on increases. The 
drop in (A^*) / (AO at critical T, would be completely counter- 
intuitive in an equilibrium picture. It would hence provide 
evidence that non-equilibrium effects such as supercool- 
ing, where such a drop would be expected, are at play. 

One difficulty of this approach is that, unlike for stable 
particles, pseudo-elastic processes such as Ntt — > N* — > 
Nn and post-decay N* — > Nn scattering of decay products 
in matter could potentially considerably alter the observ- 
able final ratio of detectable N* to N. The combined effect 
of rescattering and regeneration has not been well under- 
stood. We have argued that the formation of additional de- 
tectable resonances is negligible [2 3], while scattering of 
decay products can decrease the visible resonance yields 
except for sudden hadronization case. Transport simula- 
tions |29l confirm that rescattering generally dominates 
over regeneration in a hadronic medium. This is not sur- 
prising: If the number of reinteractions is small, rescatter- 
ing generally dominates. If the number of reinteractions is 
large, detailed balance would mean K* /K would reequi- 
librate from the higher chemical freeze-out to the lower 
thermal freeze-out temperature. 

To obtain a further quantitative estimate, we then evolved 
in time the ratios using a model which combines an aver- 
age rescattering cross-section with dilution due to a con- 
stant collective expansion. In this model, the initial reso- 
nances decay with width r through the process N* — > Nd^. 
Their decay products (No) {(No) (t = 0) = 0) then undergo 
rescattering at a rate proportional to the medium's density 
as well as the average rescattering rate. 




4 6 8 10 12 14 16 18 20 

s 1/2 [GeV] 

Fig. 3. Same as Fig. [2]more resonances 



models is robust. Namely, within the chemical equilibrium 
model the temperature of chemical freeze-out must steadily 
increase and so does the (N*) / (AO ratio. For the chemical 
non-equilibrium model the (A^*) / (AO dip primarily relies 
on the response of T to the degree of chemical equilibra- 
tion: prior to chemical equilibrium for the valance quark 
abundance, at a relatively low reaction energy, the freeze- 
out temperature T is relatively high. At a critical energy, 
T drops as the hadron yields move to or even exceed light 
quark chemical equilibrium, yet reaction energy is still not 
too large, and thus the baryon density is high and meson 
yield low. As reaction energy increases further, T increases 



d(N*) 

dt 
d(N D ) 

dt 



-r(N*), 



ri^-D^io-DjVD^Pj 



Ro + vt 



(2) 
(3) 



v is the expansion velocity, Ro is the hadronization radius, 
Pj = rij(mj, T) the initial hadron density and (cr d jVd j) is 
the hadron average flow and interaction cross-section. 

In this calculation, we have neglected the regeneration 
term R cc (cr ! ^. EL VDi)pi, since detectable regenerated reso- 
nances need to be real (close to mass-shell) particles. Fig- 
ure |4] shows how the ratios of E* I A and K* /K evolve with 
varying hadronization time within this model. It is apparent 
that measuring two such ratios simultaneously gives both 
the hadronization temperature and the time during which 
rescattering is a significant effect. Figure [4] shows the ap- 
plication of this method to K*, A(1520) and IT" H30I31L 

Assuming a long lived hadron phase, the energy depen- 
dence of most of the resonance ratios considered here has 
been calculated in a hadronic quantum molecular dynam- 
ics model. The result (figure 7 and 8 in ||29l ) is qualita- 
tively similar to the chemical equilibrium results for reso- 
nance ratios, we see a smooth rise with energy. Thus, in the 
case of chemical equilibrium, with a considerable separa- 
tion between chemical and thermal freeze-out inherent in 
Ref . [ 29 ] rescattering and regeneration will affect the quan- 
titative (N*) I (N) ratio, but will not alter the dependence 
on heavy ion reaction energy shown in figure [2] 
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Fig. 4. Two (A 7 *) / (N) diagrams over a mesh of temperature and 
lifespan. Straight lines give experimental results B30I31I . 



Thus, while at a given energy the role of post-freezeout 
rescattering/regeneration is still somewhat ambigous, an 
energy scan (experimentally conducted or planned around 
several labs [23 32 33 34 1) of resonance production should 
ameliorate some of this ambiguity. In the next section we 
will explain how fluctuations could make this ambiguity be 
experimentally resolvable at each given -\fs too. 



3 Fluctuations 

Before we discuss the significance of fluctuations in a sta- 
tistical model, we need to discuss some issues relevant 
to fluctuations and higher cumulants of particle distribu- 
tions. Generic fluctuation observables are non-trivially de- 
pendent both on non-statistical fluctuations in the systems 
volume, which can not be described by the statistical model, 
and on acceptance limitations of the detector. 

Regarding the first issue, due to our incomplete under- 
standing of how the initial state and dynamics contribute 
to fluctuations at freeze-out, the best approach is to choose 
an observable which is insensitive to any volume fluctua- 
tions. In the thermodynamic limit, where volume becomes 
a proportionality constant at the level of the partition func- 
tion, a tempting observable is the scaled variance of the ra- 
tio of two particle multiplicities measured event by event. 
That this is in fact a good guess, and all dependence of vol- 
ume factorizes, can be proven in the thermodynamic limit 
to order ~ {[AV) / (V) 2 ^ Q. The residual dependence of 
o~njn^ on the average volume (V) can be in turn elim- 
inated, in the grand canonical ensemble, by focusing on 

f = (Ni) °~dijl N2 > wnere (Ni) an d °~dyn are to be measured 
within the same acceptance. Note that this independence is 
specific to the grand-canonical ensemble, so should not ap- 
ply to scenarios where the "enhancement of strangeness" 
in A-A collisions is due to the transition between the canon- 
ical limit in elementary processes and the grand canonical 
limit in A-A H221 . In this scenario, the "strangeness corre- 
lation volume" should regulate f as in f35l . 

A different problem is the effect of a detectors lim- 
ited acceptance ( Particle (mis)identification, Limited ra- 
pidity and momentum resolution, momentum cuts neces- 
sary to eliminate jets etc) on fluctuation observables. These 
are much more difficult to model than averages, and once 
again an observable needs to be constructed insensitive to 
them. Hence, the necessity of mixed event subtraction 1 36 1 . 
Mixed events here, are defined as events where no physical 
correlation from the original event are left in. This means 
that any correlation seen is due to the imperfection of the 
detector (the fact that detector acceptance excludes parti- 
cles of pj > 1 GeV for both event A and B creates a 
correlation between A and B). To a good approximation, 
examined in the next paragraph, the measured fluctuation 

is °" 2 = ^physics + Acceptance and the mixed eVent 0ne is 

At, = trivial + Acceptance- Therefore, concentrating on 
Ayn ~ °" 2 _ Anix snou ^d eliminate acceptance effects [36 1. 

Particle abundances and fluctuations can be calculated 
from the first and second derivatives a particle's partition 
function, and related to the fluctuation of a ratio via |[T9l 

(gWg) ((^2) 2 ) (AN1AN2) ... 

«" ( Nl ) 2 + (N 2 ) 2 (Nr)(N 2 ) W 

2 _ 1 1 
0~mix ~ 777T2 /at \2 ^ ^ 

{NiY {N 2 } 

As we discussed before, provided the chemical parameters 
do not change across systems, f^', Nl = (N\ ) crj „ should 
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Fig. 6. K* /K implied from fluctuations and directly measured. 
Fluctuations and resonances are respectively from 113 8 1301 



be strictly independent of centrality and system size. This 
requirement is not satisfied by the SPS scan [39| since 
there fj./T does vary considerably. It is however satisfied by 
the RHIC upper energies. Thus, the Grand-Canonical sta- 
tistical model predicts a flat dependence with system size 
at these energies. The canonical model, on the other hand, 
would predict a "kink" within the same centrality as when 
the strangeness correlation volume approaches the thermo- 
dynamic limit 1 22 5 1. Note that a scan in centrality and sys- 
tem size within the same energy should however produce 
the same approximately horizontal bands in ¥ as are seen 
at RHIC, since thermal parameters are to a good approxi- 
mation independent of system size to the lowest centrality 
|40|: Each value of yfs, and hence p/T, should produce a 
band when scanned in centrality and system size. The re- 
sults can be seen on Fig. [5] None of the models come out 
perfectly, but the Grand-Canonical model is qualitatively 
much more similar to the data than the canonical one, as 
no kink is visible there. The discrepancy in ¥ between 200 
GeV and 62 GeV, the slight upward trend in centrality, and 
the lack of scaling between A-A and Cu-Cu should how- 
ever be closely watched, as the slight increase of fluctu- 
ations with multiplicity can not be accounted for by any 
statistical model, although one has to be careful as it could 
be an artifact of the event mixing procedure [5 1. 

Quantitatively, as reported earlier, ¥ is modeled much 
better with the inclusion of the light-quark non-equilibrium 



parameter y q , due to Bose-Einstein enhancement of fluctu- 
ations |4j . It remains to be seen weather the measurement 
of fluctuations of more particles (Fig.|7]top panel) will cor- 
roborate this conclusion. 

We now turn to a potentially important use of fluctu- 
ations, the constraint on resonance reinteraction between 
chemical freeze-out and thermal freeze-out [20 1 . The for- 
mer can be estimated, to a good approximation, by com- 
paring the fluctuations of same-charged particles (uncorre- 
lated) and opposite charged particles (correlated by K*) 

(N K -) ~ A K -"<- K *l«-> W 

Comparing the fluctuation estimate of ^ to a direct mea- 
surement should yield the amount of K°* destroyed by rescat- 
tering or regenerated through pseudo-elastic interactions. 

This analysis is shown in Fig.|6]with resonance values 
are taken from QUI . Rescattering and regeneration have,within 
error bar little or no effect on the final abundance of K*s. 
Either chemical and thermal freeze-out proceed very close 
to each other, so the amount of reinteraction is negligi- 
ble, or rescattering and regeneration of detectable reso- 
nances cancel each other out to the degree of approxima- 
tion allowed by the error bars (10 - 15%). This margin ap- 
pears somewhat below the estimates from transport mod- 
els, which are ~ 40% ||29l , and inline with fits assuming 
sudden freezeout [27 28 1. It would be very interesting to 
investigate weather a transport model [37 1 could be tuned 
to simultaneously reproduce the fluctuation inference and 
direct measurement of . 

Unfortunately, the baroqueness of the resonance decay 
tree severely limits the feasibilness of such graphic meth- 
ods, as the bottom panel of Fig [7] shows. K, n is a special 
pair of particles in that there is only one type of resonances 
that decays into both, the K*, and its lightest state is consid- 
erably lighter than the heavier ones. No similar definition 
is possible for 27* /A,p/n and (p/K, since the resonance de- 
cays equally into all pairs of decay products. For other res- 
onances, cross-contamination destroys any value of fluctu- 
ations as a graphic tool of reinteraction time. This does not 
mean, however, that such resonances are useless for con- 
straining reinteraction time, since both fluctuations [4 1 and 
resonances |2 3 1 provide a tight constraint on all statistical 
models. If a statistical model consistently describes fluc- 
tuations, but over-predicts resonance abundances, it could 
be taken as an indication of a long reinteraction times with 
detailed balance. A casual look at experimental data [ 30 1 
shows this does not seem to be the case, since most reso- 
nances are under-predictedby statistical models. 

In conclusion, we have provided an overview of the 
current results in the theoretical interpretation of both reso- 
nances and short-lived rations within the statistical model. 
We discussed how these observables should scale with dif- 
ferent energies and system sizes within various statistical 
models available on the market. We concluded that the en- 
ergy dependence of the ratios of the type (TV*) / (TV), as 
well as a simultaneus description of fluctuations and par- 
ticle yields, can be very useful for clarifying present theo- 
retical uncertainities. 
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Fig. 7. Equilibrium and non-equilibrium predictions from fluctu- 
ations (top panel) and resonances (bottom panel) with the statis- 
tical parameters taken from [24|. In the bottom panel, Full sym- 
bols show the (N*) IN ratio inferred from the correlation in the 
two models, while empty symbols show, wherever possible, the 
estimate from a dlJ „ comparisons 
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